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ABSTRACT 

It is attractive to suppose for several astrophysical reasons that the universe has cfose 
to the critical density in light (~ 30 eV) neutrinos which decay radiatively with a 
lifetime of ~ 10^^ sec. In such a cosmology the universe is reionized early and the 
last scattering surface of the cosmic microwave background significantly broadened. 
We calculate the resulting angular power spectrum of temperature fluctuations in 
the cosmic microwave background. As expected the acoustic peaks are significantly 
damped relative to the standard case. This would allow a definitive test of the decaying 
neutrino cosmology with the forthcoming MAP and PLANCK surveyor missions. 
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1 INTRODUCTION 



al. 1998; Garnavich et al. 1998| ) actually suggest Q, 



The absence of Gunn-Peterso n absorption by diffuse n eutral 
hydrogen in quasar spectra ( Gunn fc Peterson 1965| ) indi- 
cates that the intergalactic med ium (IGM) is highly io nized, 
out to a redshift of at least 5 (Giallongo et al. 1994). The 
integrated flux of ionizing UV photons due to quasar s them- 
selves is believed to be inadequ ate for this purpose (Giroux 
& Shapiro 1996; Madau 1998). It has been proposed that 



the required flux can arise from the radiative decays of relic 
neutrinos of mass ~ 30 eV which would also provide t he crit- 
ical d ensity in dark matter ( Rephaeli fc Szalay 1981 ; Sciama 
|l982 ). This hypothesis can naturally acc ount for the tem- 
perature of Lyman-Qf clouds (Sciama 1991) as well as explain 
an anomaly found in the abundance of He I in three high- 
r edshift Lyma n-limit systems of the QSO HS 1700-f6416 
(sciama 1994). Moreover, the decays of the neutrinos in 
galactic halos such as ours can account for the observed 



0.3-0.4 

(with a possible cosmological constant). The decaying neu- 
trino cosmology can, in principle, accomodate this as long 
the constraint (|l]) is respected, i.e. for a high value of h.) Di- 
rect searches for the expected UV line at — 13.7±0.1 eV 
(i.e. A = 905 ±7 "A) have so far proved unsuccessful but new 
results are expected soon from the EURD de tector on the 
Spanish MINIS AT 01 satellite ( ^ciama 1997b| ). 

As is well known, large-scale structure formation is 
problematic in a HDM cosmology because the primordial 
density perturbation is seve rely damped on small scales du e 
to neutrino free-streaming (Bond, Efstathiou & Silk 1980). 
This would lead to supercluster-size objects forming first 
and smaller objects forming subsequently through their frag- 
mentation, wh ich is in conflict with a variety of observations 
(Peebles 1993). However the problem may be solved if there 



is an additional source of small-scale fluctuations, e.g. from 



ionizati pn of the interstellar medium, which again is di ffi 



cult to l account for by conventional oourcoo ( jciama 1003 )^ 
Variou a aotrophyoical and coomological conotrainto can then 
down the pBrameters nf the Hec Ryinp; hnt riBrk matter 



a network of cosmic topological defects ( Villumsen, Scherrer 
& Bertschinger 1991; [Gratsias et al. 1993|) 



pm 



It had been noted some time ago (Bond fc Efstathiou 



(HDM) cosmology rather precisely (Sciama 1997a) 



= 27.4±0.2eV ^ n^ft= 0.293 ± 0.003, (1) 
~ (1 - 2) X 10^^ sec, (2) 
h = 0.548 ± 0.003, for Q„ = 1. 

(Here n^, and denote the fraction of the critical density 



1984) that reionization of the IGM, e.g. by early star forma- 
tion, leads to suppression of the anisotropy in the cosmic mi- 
crowave background (CMB) on angular scales smaller than 
the horizon size at the reionization epoch. This would then 
appear to provide a potential test of the decaying HDM cos- 
mology. Although much work has been done on reionization 
in cold dark matter (CD M) and (isocurvature) baryonic dark 
matter (BDM) mo d els (|Sugiyama, Silk fc Vittorio 1993^ ^ 



ble con itant in units of 100 km sec Mpc . We note that 



in neut rinos and baryons, respectively, and h is the Hub- fc Sugiyam a 1994; Dodelson fc Jubas 19951; 



De Bernardis 



most dynamical esti mates of the c lustered mass density upto 
supercluster scales ( Peebles 199S ) as well as r ecent observa- 
tions of Type la supernovae at high redshift ( Perlmutter et 



et al. 1997 ), these results are not directly applicable to the 
present case where the universe is gradually reionized by the 



decaying neutrinos. (An exception is the work of Dodelson 
fc Jubas (1995) who considered a variety of possible ion- 
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ization histories and estimated the CMB fluctuation signal 
expect ed in speciflc experiments.) Hence it is necessary to 



make a specific study of the problem. Such a calculation 



has been performed already by Tuluie, Matzner and Anni- 
nos (1995) (see also Anninos et al. 1991) who numerically 
simulated the growth of structure in a HDM universe and 
followed photon trajectories in the reionized universe in or- 
der to construct CMB temperature maps. However these au- 
thors fixed the amplitude of the initial density perturbation 
by requiring that the autocorrelation function of the dark 
matter particles should equal the observed value for clusters 
(C(^) = (^/''o)~^'*) at a redshift of 2 = for separations 
r ~ 25/1"^ Mpc. Since the COBE discovery of super-horizon 
scale CMB fiuctuations this is no longer appropriate and 
the primordial density perturbation should be normalized 
via the COBE quadrupole. It is then seen that the HDM 
model has too little power on cluster and smaller scales to 
match observations and, as mentioned earlier, some other 
source of small-scale power must be sought. Since this issue 
is model-dependent, the effects of reionization on the CMB 
anisotropy ought to be computed without speciflc reference 
to the manner in which structure was generated. 

In this Letter we flrst make the necessary connection 
between the ionization history of the IGM and the lifetime 
of the decaying neutrinos. We then incorporate this in a nu- 
merical 



licode which compuLes Lhe CMB auisoLiopy with high 
precisi an. To compare with observational data, we normalize 
the primordial density perturbation to COBE and present 
the results as an angular power spectrum. 



2 REIONIZATION IN THE DECAYING 
NEUTRINO COSMOLOGY 

In the standard cosm ology the universe (re)combines at a 
redshift of Zrec ~ 1300 ( [Peebles 1968, 198l| ). In particular the 
last scattering surface (LSS) of the relic blackbody photons 
is well approximated by a gaussian of wi dth A z — 78 located 



at a redshift z ~ 1065 (Jones & Wyse 1981:). The presence 



of decaying neutrinos generates UV photons which reion- 



ize the IGM, thus si gniflcantly broadening the LSS (scott 



Rees & Sciama 1991). We study the ionization history from 
early times by direct integration of the rate equation for the 
ionization fr action of hydrogen, x h = n^/nu. following pre- 



vious work (Salati fc Wallet 1984 ; Asselin, Girardi & Salati 
198^; podelson fc Jubas 19921). 



At high temperatures both helium and hydrogen are 
fully ionized and xh ~ 1. Helium (re)combines while hy- 
drogen is still fully ionized and this process can be approx- 
imated by the Saha ionization equilibrium equation, quite 
independently of the subsequent evolution. The evolution of 



the free electron fraction through hydrogen (re)combination 



and rei( )nization is driven by two distinct processes — tran- 



sitions via excited states, and transitions directly to the 
ground state. 



da:H 
dt 



rig nisUi 

= a I ais (Ti 

riH riH 



(3) 



+ aC [ a— - »^^-Bi/kTB 
nil 



„nis _ 
■ p e 



where Q13 and a axe the recombination rates to the Is level 



and excited states respectively, (3 is the ioniza tion rate from 
theground state, C is a correction factor (Peebles 1968 



1981) and Bi — 13.6 eV is the ground state binding en- 
ergy. In the standard picture of (re)combination of hydrogen, 
the first component is negligible because every combining 
electron-proton pair releases a photon which immediately 
ionizes another neutral hydrogen atom. However it cannot 
be neglected when there is an an additional injection of ion- 
izing photons from the decaying neutrinos at a rate 

/dni\ _ 

V di / decay 

where ni and are the number density of ionizing photons 
and neutrinos respectively and is the neutrino lifetime. 
Then, including the effect of dilution due to expansion and 
the redshifting below the Rydberg energy, ni evolves in time 
according to 



(4) 



dni 



+ 



2 , 

aiaUg -\ 



SO the equilibrium value of ni is 



(ni)c 



[aims + {3 + Bi/kT)HY 



(5) 



(6) 



Except for a small region when sh is clos e to unity, ni is well 
appr oximated by its equilibrium value (Dodelson & Jubag 



1992). It is also a good approximation to replace the number 
density nis of hydrogen atoms in the Is state by (1 — a::H)wH. 
Equation Xfj then reads 



da:H 
' dt 



(7) 



aisXiinn - (Ti(l - xn) 
+ aC (dXnnH — /?(1 — a;H)e 



Gi(\ - XH)nH + {3 + Bi/kT)H 

Bi/kTB\ 



In addition we must follow the evolution of the matter and 
radiation temperatures. We modify the standard ev olution 
as incorporated in the computer code COSMICS (Ma & 
Bcrtschinger 1995) by including the energy, E^, — Bi, re- 



leased by the decaying neutrinos (as long as there is neutral 
hydrogen present to capture the photons). 

We integrate Equation (M) numerically using a 
semi-implicit method. The resulting ionization fraction, 
optical depth and cumulative visibility function (= 
e~'^ (dt/d z) dz) are shown in Figure ^ . For the Hub- 
ble parameter we use h — 0.55 (Sciama 1997c). The baryon 
density obtains from taking D.Bh'' = 7 x 10"'^ as indicated by 
the primordial abundance of ^He inferred from observations 
of blue compact galaxies by Olive, Skillman & Steigman 
(1997), the h igh D abundance measured in quasar absorp- 
tion systems ( Bongaila, Wampler & Cowie 1997 ; Webb et al. 
1997 ) and the abundance of ' Li in Pop II stars ( ponifacio 
& Molaro 1997). If however we adopt the higher *He abun- 



dance obtained by Izotov, Thuan & Lipovetsky (1998) in 
conjunction with the lower D abundance of Buries & Tytler 
(1998) then the appropriate value is Qb/i^ = 1.88 x 10~^. 
Both values for the baryon density, viz. nB = 0.023, 0.062 
yield a good for th e fit to the theoretical expectations 
( Fiorcntini et al. 1998| ). We also consider two values for the 
neutrino lifetime, t^, — 10^^, 2 x 10^^ sec. It is seen that 
the IGM is fully ionized for a redshift below ~ 30 although 
the optical depth reaches unity only at a redshift exceed- 
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Figure 1. The fractional ionization, optical depth and cumula- 
tive visibility function in the decaying neutrino cosmology for a 
neutrino lifetime of (a) 10^^ sec and (b) 2 X 10^^ sec. 
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Figure 2. Angular power spectrum of CMB anisotropy in the 
decaying neutrino cosmology for (a) = 0.023 and (b) = 
0.062, with T„ = 10^^, 2 X lO^^ sec. 



ing ~ 1000. The LSS is significantly broadened as is evident 
from the slow growth of t he cumulative visibility func tion. 
It had been argued (Zeldovich & Sunyaev 196S;) that 



the universe must have been neutral during the epoch 



soon 



-7/9 



< 2 < 2rcc in order not to induce a "y" distor- 



tion in the blackbody spectrum of the CMB. However this 
assumed that all matter is constituted of baryons. Allowing 
for non-baryonic dark matter, this argument is ev aded and 
the uni verse is n ot required to have (re)combined (Bartlett 
& Steb jins 1991). We have checked that there is no conflict 



in the decaying HDM cosmology with even the very restric- 
tive bound ^< 2.5 X 10"^ from the CORE FIRAS data 



(Fixsen et al. 1991) 



3 CALCULATION OF THE CMB 
ANISOTROPY 

As mentioned earlier, the anisotropy in the CMB will be 
damped as a result of the broadening of the last sca ttering 
surface. To calculate this we use the CMBFAST code (Seljak 



& Zaldarriaga 1996) incorporating the ionization fraction 



evolution described above. The primordial density pertur- 
bation is assumed to be scale-invariant and is normalized on 
super-horizon scales to the CMB qu a.drupole moment me a- 
sured by COBE: Qrms-ps ^ 18 fiK ( [Bennett et al. 1996| ). 

The calculated angular power spectrum is shown in Fig- 
ure ^. As expected the acoustic peaks are damped by reion- 
ization, with the extent of damping decreasing as the neu- 



trino lifetime increases. (We stress that the damping effect 
is largely due, not to the reionization itself, but rather the 
presence of more free electrons at all reshifts which widens 
the last scattering surface. In particular, sudden reionization 
at a late redshift would not result in significant damping.) 
An increase in the assumed baryon density alters the rela- 
tive level of damping of the second and third acoustic peaks. 
A selection of data from recent ground-based and balloon- 
borne experiments is also shown for comparison. 

To obtain a quantitative measure of the (dis)agreement 
of the decaying neutrino predi ction with present ob serva- 
tions we perform a analysis (Lineweaver et al. 1997) after 
convolving the predicted CMB power spectrum with the ex- 
perimental window functions. The for the decaying HDM 
model ranges between 500 — 700 for the allowed range of 
Ti, and Qb, for 29 degrees of freedom (32 data points - 2 fit- 
ted parameters - 1 COBE normalization) . Since there is less 
than 0.5% probability of obtaining a exceeding 52.34, the 
decaying HDM theory clearly does not fit the trend in the 
present data. However it would be premature to draw any 
strong conclusion from this, given that there may still be 
large systematic uncertainties in the present data. A defini- 
tive test will however be possible with the forthcoming MAP 
and PLANCK all-sky surveyor satellite missions, which are 
expected to reliably measure t he small-scal e anisotropy with 
an accuracy of a few per cent ( gmoot 1997 ). Moreover these 
experiments will also measure the CMB polarization, thus 
provi ding an independeri t constraint on the ionization his- 
tory dKcating et al. 1998|) . 
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4 CONCLUSIONS 

We have computed, in a consistent manner, the angular 
power spectrum of CMB anisotropy to be expected in the de- 
caying HDM cosmology. The damping of the acoustic peaks 
due to the gradual reionization of the IGM is shown to pro- 
vide a clear test of the model. 

Other possible observational signatures of reionization 
have been considered. With regard to the CMB, although 
large angular-scale fluctuations are damped, new fluctua- 
tions are generated on arcminute-scal es by the motions of 
the i o nized gas — th e Vishniac effect ( Ostriker fc Vishniac 
198(: ; Vishniac 1987 ) . This has been studied in some detail 



for both CDM and EDM models and various constraints in- 
ferred fro m observational l imits on small angular - scale fluc- 
tuations dsfstathiou IQSSl IHu, Scott fc Silk 1994 jPodelson 
fc Jubas 199£| ; |Hu fc White 1996] ). For HDM, the small- 



scale power in the primordial (inflationary) density pertur- 
bation is generically suppressed which would normally imply 
a negligible Vishniac effect. However as we have noted some 
additional source of small-scale power such as cosmic topo- 
logical defects is needed in any case to allow galaxies to form 
in an HDM universe. Lacking a precise knowledge of such 
fluctuations, it is not possible to make a model-independent 
calculation of the Vishniac effect in the present case. 



ACKNOWLEDGMENTS 

We thank Charlie Lineweaver for providing the CMB data 
and experimental window functions, Ed Bertschinger for the 
COSMICS code and Uros Seljak and Matias Zaldarriaga for 
the CMBFAST code. D.W.S. is grateful to M.U.R.S.T. for 
their financial support of this work. 



Hu W., White M., 1996, A&A, 315, 33 

Jones B.J.T, Wyse, R.F.G., A&A, 1985, 149, 144 

Izotov Y.I., Thuan T.X., Lipovetsky V.A., 1998, ApJ, 500, 188 

Keating B., Timbie P., Polnarev A., Steinberger J., ApJ 495, 580 

Lineweaver, C.H., Barbosa, D., Blanchard, A., Bartlett, J.G. 

A&A, 1997, 322, 365L 
Ma C.-P., Bertschinger E.. 1 995. ApJ, 455, 7 
Madau P, a,stro-ph/9804280 

Olive K.A., Skillman E., Steigman G., 1997, ApJ, 483, 788 
Ostriker J. P., Vishniac E.T., 1986, ApJ, 306, L51 
Peebles P.J.E., 1968, ApJ, 153, 1; 1981, ApJ, 248, 885. 
Peebles P.J.E., 1993, Principles of Physical Cosmology, Princeton 

University Press 
Perlmutter S. et al., Nat, 391, 51 

Rephaeli Y. Szalay A.S., 1981, Phys. Lett. B, 106, 73 
Salati P., Wallet J.C., 1984, Phys. Lett. B, 144, 61 
Sciama D.W., 1982, MNRAS, 198, Ip 
Sciama D.W., 1991, ApJ, 367, L39 

Sciama D.W., 1993, Modern Cosmology and the Dark Matter 

Problem (Cambridge University Press) 
Sciama D.W., 1994, ApJ, 422, L49; 426, L65 
Sciama D.W., 1997a, MNRAS. 289. 945 
Sciama D.W., 1997b, |astro-ph/9709243| 
Scott D., Rees M.J., Sciama D.W., 1991, A & A, 250, 295 
Seljak, U., Zaldarr iaga. M.. 1996. Ap J. 469, 437 
Smoot G.F, 1997, [istro-ph/9705135 

Songaila A., Wamplcr E. J., Cowie L.L., 1997, Nat, 385, 137 
Sugiyama N., Silk J., Vittorio N., 1993, ApJ, 419, LI 
Tuluie, R., Matzner, R.A., Anninos, P., 1995, ApJ, 446, 447 
Villumsen J.V., Scherrer R., Bertschinger E., 1991, ApJ, 367, 37 

(Erratum, 381, 601) 
Vishniac E.T., 1987, ApJ, 322, 597 
Webb J K et al., 1997, Nat, 388, 250 
Zeldovich Ya B., Sunyaev R.A., 1971, Ap&SS 4, 301 



REFERENCES 

Anninos P., Matznor R.A., Tuluie R., Centrella J., 1991, ApJ, 
382, 71 

Asselin X., Girardi G., Salati P., 1988, Nucl. Phys. B, 310, 669 

Bartlett J.G., Stebbins A., 1991, ApJ, 371, 8 

Bennett C.L. et al., 1996, ApJ, 464, LI 

Bond J.R., Efstathiou G., 1984, ApJ 285, L45 

Bond J.R., Efstathiou G., 1987, MNRAS 226, 655 

Bond J.R., Efstathiou G., 1980, Phys. Rev. Lett. 45, 1980 

Bonifacio P., Molaro P., 1997. MNRAS. 285. 847 



Buries S., Tytler D., 1998, s,stro-ph/9803071 



De Bernadis P., Balbi A., De Gasperis, G., Melchiorri A., Vittorio 
N., 1997, ApJ, 480, 1 

Dodelson S., Jubas J.M., 1992, Phys. Rev. D, 45, 1076 

Dodelson S., Jubas J.M., 1995, ApJ, 439, 503 

Efstathiou G., 1988, Large Scale Motions in the Universe: A Vat- 
ican Study Week, Rubin V.C. & Coyne G.V., eds., Princeton 
University Press, Prin ceton, n.299 



Fiorentini G. et al., 1998, |astro-ph/9803177 
Fixsen D.J. et al., 1996, ApJ, 473, 576 
Garnavich, P.M. et al., ApJ, 493, L53 
Giallongo E. et al., 1994, ApJ, 425, LI 
Giroux M.L., Shapiro P.R., 1996, ApJS, 102, 191 
Gratsias J., Scherrer R., Steigman G., Villumsen J.V., 1993, ApJ, 
405, 30 

Gunn J.E., Peterson B.A., 1965, ApJ, 142, 1633 

Hu W., Scott D., Silk J., 1994, Phys. Rev. D, 49, 648 

Hu W., Sugiyama N., 1994, ApJ, 436, 456 



© 1998 RAS, MNRAS 000, 000-000 



